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ABSTRACT 

From the 389 2002 OGLE-III observations of Galactic Bulge microlensing events we 
select 321 that are well described by a point-source point-lens lightcurve model. From 
this sample we identify n = 1 event, 2002-BLG-055, which we regard as a strong plan- 
etary lensing candidate, and another, 2002-BLG-140, which is a possible candidate. If 
each of the 321 lens stars has 1 planet with a mass ratio q — m/M = 10~ 3 and orbit 
radius a = iig, the Einstein ring radius, analysis of detection efficiencies indicates that 
14 planets should have been detectable with A% 2 > 25. Assuming our candidate is 
due to planetary lensing, then the abundance of planets with q = 10~ 3 and a = Re is 
rip « n/14 — 7%. Conversion to physical units (Mj up , and AU) gives the abundance 
of 'cool Jupiters' [m « Mj up ,a « 4 AU) per lens star as n p w n/5.5 = 18%. The 
detection probability scales roughly with q and (Ax 2 ) -1 / 2 , and drops off from a peak 
at a w 4 AU like a Gaussian with a dispersion of 0.4 dex. 
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1 INTRODUCTION 
1.1 Microlensing 

Gravitational lensing was first proposed by [Einstein! jl936f) 
- a consequence of his theory of general relativity, it is the 
effect of a massive body warping space-time and bending 
light around it, acting as a lens. This produces two or more 
images of background objects distorted around the 'Einstein 
ring', which can be thought of as a scale radius (Re) of the 
lens and depends on its mass, M. The angular radius of the 
ring is given by: 



Re 
Dl 



D s - D L AGM 



D s Dl 



(1) 



where Dl and Ds are the observer - lens and observer - 
source distances respectively. For objects of stellar mass, 
this radius is below the resolution limit of our telescopes, 
and we observe the phenomenon of microlensing. In this 
case, we cannot resolve multiple images of the source, but 
detect an increase in brightness as the individual amplifi- 
cations combine to produce a brighter stellar image. As we 
observe a foreground (lens) star passing in front of a back- 
ground (source) star, we see an increase in brightness up 
to a maximum amplitude Aq at the time to of closest ap- 
proach, followed by a decrease. This gives a symmetrical 



light curve around to, providing the timescale for the event 
is short enough for parallax effects to be negligible. The 
timescale depends on the lens transverse velocity, relative 
to the observer-source line of sight, v, and is given by the 
length of time it takes the source to cross the Einstein ring 
of the lens: 

te = ^. (2) 

v 

This depends on both mass and distances, through the 
dependence on Re- v is also indirectly influenced by dis- 
tance, as the positions within the galaxy of the source and 
lens stars will give them average motions based on the or- 
bital velocity at their galactic radius. For a survey looking 
toward the galactic bulge (i.e. (Ds — Dl) / DsDl ~ kpc), 
a microlensing event due to a stellar mass (M = Mg) lens 
will last a few weeks. A Jupiter mass object will have tE ~ 
1 day, while an Earth would produce an effect which would 
be observable for only a few hours. 

If a stellar lens has a planetary companion, this will 
appear as a brief deviation from the stellar microlens- 
ing light curve, wh ere the am plificati ons due to the star 
and planet combine <Mao fc Paczvriskil99lHGould fc Loebl 
Il992t [Bennett fc Rhieil99fJ) . These deviations should be ob- 
servable, even for low mass objects, and therefore provide 
a method to detect extra-solar planets which is potentially 



2 Snodgrass et al. 



more sensitive to small planets than radial velocity searches, 
given a programme of intensive monitoring with data taken 
many times per hour. 



1.2 The OGLE project 

The main difficulty with using gravitational microlensing 
to search for exoplanets is the low probability that any 
given source star will be lensed. This can be overcome by a 
monitoring programme which studies dense star fields and 
can detect a flux increase in any of the stars in the field. 
The Optical Gravit ational Lensing Experiment (OGLE) 
JUdalski et alJll994) monitors ~ f50 million stars in the 
galactic bulge, and alerts the astronomical community via 
Internet and e-mail alerts as soon as a lensing event is seen 
to begin. 

OGLE began in its current form, OGLE-III, in 2001. 
This version uses Dif ference Image Analysis (DIA) photom- 
etry dWozniakllioOol) . instead of the profile fitting (PSF) 
method used for the previous three seasons. This has sub- 
stantially increased the detection rate of stellar microlensing 
events by OGLE - there were 167 events discovered during 
the three years of OGLE-II, compared to 389 recorded dur- 
ing the 2002 season with OGLE-III. 

We use OGLE 2002 data to put firm upper limits on the 
population of planets in the bulge o f our galaxy. Our method 
follows that of lTsapras et alj ll2003f) . and is described in Sec- 
tion|21 before Section|3]looks at the results found for the 2002 
data. The abundance of planets is discussed in Section [I] 



2 LIGHT CURVE FITTING 
2.1 PSPL Model 

A Point-Source Point-Lens (PSPL) model of microlensing 
can describe the light curve entirely based on four param- 
eters; to and £e as defined in Section [l.ll the baseline (un- 
lensed source) magnitude Jo, and the impact parameter uo- 
The angular separation of the unlensed source and lens, mea- 
sured in units of 6e, varies in time, 



u(t) 
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(3) 



and relates to the amplification of the signal, which governs 
the shape of the light curve with (Fig. 12. 1 1 time, by: 



A(u(t)) 



u + 2 

Uyju 2 + 4 



(4) 



The impact parameter ito is the separation of the source 
and lens at closest approach, which is at to, and gives the 
maximum amplification Ao = A(uo). 

PLENS is a program which fits these parameters to mi- 
crolensing data by minimising the \ 2 measure of 'badness 
of fit'. In the case of microlensing, where h and ai are the 
reported magnitude and error bar for each data point, 
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Figure 1. Light curve and normalised residuals for event 2002- 
BLG-024 - an example of a good PSPL fit. 



J(t) = J - 2.5 log (A(t) + b) 



(6) 



where the source magnification A(t) is defined in Eqn.^] Jo 
is the source star's baseline magnitude, and b accounts for 
the additional flux from other stars that are blended with 
the source. 

This method produces reasonable fits to the microlens- 
ing peak, however it was often found that the residuals in the 
flat (baseline) parts of the curve were larger than the pho- 
tometric errors <Ji reported for data points in this region. 
This was overcome by introducing a further fit parameter, 
k„, which scales the error bars. However this scaling often 
causes the error bars around the central peak to become 
much larger than the residuals, and would mask any inter- 
esting discrepancies. 

The scatter around the baseline can be explained by the 
fact that the fields observed are very densely populated, and 
so one further parameter was introduced: a , the crowded 
field error. When ao and fc CT are used in conjunction they 
adjust the error bars to provide a much more accurate rep- 
resentation of the residuals. For ao and at expressed in mag- 
nitudes, the adjusted error bars are then given by: 



A*(t) 



(7) 



and the best fit is found by maximising the likelihood crite- 
rion L, and is equivalent to minimising 



2 In J = x + 2 ^ In S j 



(8) 



J(t) represents the model light curve, and is given by: 



There is one more fit parameter which PLENS can ad- 
just, the blending fraction b. This describes the amount of 
the baseline flux that is due to light from non-lensed stars 
in the crowded fields studied. 



2.2 Planet detection 

Ax 2 maps (Fig. 12.21 were produced by adding a theoretical 
planet, with mass ratio q = m/M, at each point on a grid 
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blg-024 q=0.001 ky 2 =25 best Ax 2 =6.3 




covering the lens plane. The resultant x 2 value is calculated, 
using a binary lens model, for each point. This results in a 
plot of 'detection zones' around each data point, where the 
addition of a planet to the model gives |A^ 2 | > Axt, where 
A\ 2 is the difference in \ 2 between the binary and PSPL 
models, and Axt is a threshold value of A% 2 - The plots show 
a white zone where the \ 2 value is significantly increased, 
compared to a PSPL fit, by the addition of a planet. This 
means that the presence of a planet in these zones can be 
ruled out - if there were a planet at that point, it would 
have significantly affected the data point at the centre of the 
zone. A black spot on the detection map would correspond 
to A\ 2 < —Axt, implying that the planet model improves 
the fit from the PSPL by a significant margin, which would 
be strong evidence for a planet. 

The probability of detecting a planet (Fig. [3J with a 
given mass ratio, q, and circular orbit of radius a, can be 
found by: 

P(det\a,q) = J P(det\x,y,q)P(x,y\a)dxdy. (9) 

This equation simply sums, over all possible grid positions, 
the product of probability of detection at that position: 

••<**• ».«>-{; x£i: &xi <«» 

and a geometrical term giving the probability that a ran- 
domly orientated orbit, of radius a, will put the planet at 
that point: 

} if r — \/ x 2 + y 2 < a 

2xaV« 2 -r 2 V (11) 

otherwise. 

The first term comes from the A\ 2 map - it decides whether 
or not a planet at that point falls into a detection zone. 
As the detection zones are, in general, clustered around the 
Einstein ring (see Fig. 12. 21 . the highest detection probability 
is found at around a w Re- This is not surprising, as it is 
when a planetary lens interacts with a stellar lens at Re 
that its effect is greatest. 

2.3 Choice of detection threshold 

The total number of data points in the 321 OGLE 2002 
events studied is 30582. A statistical treatment of this sam- 
ple gives us the largest error expected from this much data 
(as a multiple of a) using: 
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Figure 3. Detection probability for event 2002-BLG-024. The 
curves show the probability of detection at A\ 2 thresholds of 
Ax| = 25,60 and 100. 



S F (N) « 2. 14V log TV, (12) 

assuming a Gaussian distribution of errors. For the 2002 
data, with N = 30582, the expected value of the largest 
residual is 4.53<r, corresponding to a A\ 2 ~ S% ~ 20. This 
implies that the choice of Axt = 25, corresponding to de- 
viations from the PSPL fit of more than 5<r, is appropriate 
for this data set - any signal above this threshold should be 
larger than noise, and due to a real anomaly. 

2.4 Analysis of photometric scatter 

To judge the validity of the stated (and adjusted) error bars, 
the scatter in photometric data was investigated. There ex- 
ists within the data some 14 nights of observation when more 
than 10 individual photometric measurements of the same 
source were acquired in a night's observation (e.g. 2002- 
BLG-269 - Fig. |1J. For events where there is a night of 
intensive observation away from the peak, the light curve 
would not be expected to change much over the course of 
the night. Therefore the scatter in the data points taken on 
these nights can be used to judge errors. 

Data were taken from the 128 nights, in 22 events, where 
10 or more frames were observed in a single night. The pho- 
tometric scatter, Si, was determined for each by finding the 
root-mean-squared deviation from the mean x of the mag- 
nitude data x: 

N 

S ? = J^iJ2fa-x) 2 > ( 13 ) 

3=1 

where there are N data points on the night i. This RMS 
scatter is a reasonable way to judge the errors in the pho- 
tometry, as PLENS assumes a Gaussian distribution of er- 
rors. Si, <Ji and Si correlate very well; the scatter observed 
is consistent (Fig. |KJ with the mean reported error bar for 
the nights studied - scaling directly by a factor of ~1.2. A 
histogram of the values of k„, the error bar scaling factor, 
(Fig. |SJ fitted to the 321 analysed events, agrees with this 
hypothesis - it peaks around 1.2. 



4 Snodgrass et al. 



J 


t : 

~f-5*»" -a,,, .T., J ,< 8 #\«»,|lt> 









Figure 4. 2002-BLG-269 - an example of an event with well 
sampled nights, which are clearly seen as the seven vertical lines 
of data at ^constant t (t on x/R^) in the residuals plot. They can 
be seen less clearly in the light curve plot. 



The sample of 22 events that contain these well-sampled 
nights represents 7 percent of the 2002 events studied. They 
are spread throughout the data; it would be a reasonable as- 
sumption to extend the conclusions of this scatter analysis 
to all events. This would indicate that the scatter in data in 
other events is approximately the same as the quoted pho- 
tometric error. Si and Si also correlate strongly - with the 
adjusted error bars consistently larger. We concluded that 
the scaling of error bars by PLENS is sufficient to encom- 
pass any likely photometric scatter, which should prevent 
any noise appearing as a planetary candidate. 



3 2002 OGLE EVENTS 

3.1 Identification of suitable events 

We selected events from the OGLE data that were reason- 
able approximations to a PSPL model. As the OGLE EWS 1 
system looks for any brightening which could be a lensing 
event, it occasionally puts out a false alert - the 'event' later 
turns out to be a variable star, for example. Also, a number 
of the events observed have high mass ratio binary lenses 
(i.e. a binary star acting as a lens), which PLENS is not 
designed to deal with. The non-PSPL events were removed 
from the sample. 

These events were identified by assessing the accuracy of 
the fitted parameters for each event. A Monte-Carlo method 
was used to produce data sets based on the OGLE data, with 
new values for each point chosen within a Gaussian distri- 
bution based on the stated error bars. The RMS variation 
in the parameters fitted to these data sets then provided 
an uncertainty measurement on each parameter. Fits with- 
out blending were carried out on 200 fake data sets for each 
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Figure 5. RMS Scatter (Si — crosses), reported errors (<Tj - cir- 
cles), and adjusted error bars (sj - triangles), compared to the 
mean magnitudes of the detections taken that night, I. The scat- 
ter clearly increases for fainter sources, as expected, due to the 
lower signal-to-noise ratio at lower brightness. The scatter and 
reported error bars are largely consistent, with scatter generally 
greater than the reported error bars, but within the scaled error 
bars. 



event, providing error bars with a 1/V200 w 7 percent un- 
certainty. 

It was found that the baseline magnitude Io and peak 
time to were well determined for most events, while the am- 
plification Ao and event timescale £e can vary considerably 
in events without a clearly defined microlensing light curve. 
The uncertainties in these parameters, cja q and at E , were 
used to provide filters to remove poor events. A sieve of four 
filters was used to catch non-PSPL events: 

(i) The three events listed as mistakes on the OGLE web 
page were removed. 

(ii) Events with oa > or a tE > itE were removed 
- i.e. events with over 50 percent uncertainty in either of 
these parameters. 

(iii) If an event's ga and <rt E were greater than 3 times 
their respective median values, the event was removed. This 
caught events that did not contain enough data to constrain 
these parameters, for example those only observed during 
the rising part of the light curve. 

(iv) A brief visual inspection caught those events with 
clear non-PSPL light curves, that had passed through the 
other filters. 

From the original list of 389 events the filters rejected 
3, 42, 10 and 13 events respectively - a total of 68, which 
are listed in Appendix|E| This left 321 events that were then 
searched for planetary candidates. 
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Figure 6. Histogram of values of k a , the error bar scaling factor. 
The dashed line shows the cumulative histogram. It can be seen 
that the value for is most commonly found in the range 1.1 — 
1.3, which implies that the reported error bars are smaller than 
the actual error, according to these 7-parameter fits. The median 
value of k a found here is in excellent agreement with the scaling 
factor found in Section l2.4l which implied that actual errors where 
~ 1.2 times the reported errors. 



Table 1. Data points with Ax > 25 following a fit including 
blending with q = 10 -3 . The reason for each anomaly is listed in 
the last column; those marked with an asterisk were suggested by 
Udalski (private communication). Those events listed as plane- 
tary candidates are the events for which the observation in ques- 
tion was OK, and no alternative explanation was given for the 
anomaly. 



Event 




Frame 

(HJD-2450000) 


Ax 2 


Anomaly 


2002 


BLG- 


1 a n 
14U 


z4o4.oliOo 


1 in nQ 
11U. Uo 


Planet candidate? 


2002 


BLG- 


055 


Z424.oyoD0 


64.97 


Planet candidate. 


2002 


BLG- 


205 


2472.50459 


48.06 


Clouds.* 


2002- 


BLG- 


110 


2472.50459 


47.53 


Clouds.* 


2002 


BLG- 


085 


2426.74382 


43.63 


Clouds.* 


2002 


BLG- 


321 


2517.62203 


34.60 


Lens is a stellar mass 












binary. * 


2002 


BLG- 


009 


2143.52890 


32.29 


Lone high data 












point well before 












peak. Data taken 












in 2001 observing 












season, early in 












OGLE-III when 












photometry was not 












as accurate. 


2002- 


BLG- 


186 


2472.65734 


32.01 


Clouds.* 


2002 


BLG- 


011 


2410.74028 


28.54 


Clouds.* 


2002 


BLG- 


311 


2472.65734 


28.50 


Clouds.* 


2002 


BLG- 


156 


2471.54596 


28.15 


Bad image - highly 












out of focus.* 


2002- 


BLG- 


065 


2453.54649 


25.44 


Clouds.* 



3.2 Planetary candidates 

A web page interface simplified the search for planetary can- 
didates. By listing the events in order of decreasing A\ 2 
- with Ax 2 now representing the 'best' value from each 
map; the magnitude of Ax 2 where it is at its most nega- 
tive, and the planet model is most likely to be correct - the 
events with the largest discrepancy from PSPL could be eas- 
ily identified. Examination of the original photometric data, 
the Ax 2 maps, and PLENS's output of the x-y grid position 
of the best point, allowed the data point corresponding to 
the 'planet' to be identified. 

Of the analysed sample of 321 events, the 12 events 
listed in TableQhave Ax 2 > 25. These high Ax 2 events do 
not necessarily represent planetary candidates. Examination 
of the individual frames which gave the planet-like anomalies 
(Udalski, private communication) revealed that 10 of these 
events have anomalies due to other factors, such as clouds 
etc., as detailed in Tabled 

Of the two remaining events, we regard 2002-BLG-055 
as a good planetary microlensing candidate, and 2002-BLG- 
140 as a possible candidate. 



3.2.1 2002-BLG-055 

One very high data point shortly after the peak of this well 
sampled event appears to be due to a binary lens with q 
in the expected range for a star and massive planet system 
(Fig. 0. This point was also .judged to be a g ood planetary 
candidate bv Ijaroszvnski fc Paczv hski (2002J). They find a 



Table 2. Comparison of fitted parameters to event 2002- 
BLG-055 found by our PSPL model, to those found by 
Jaroszyhski and Paczyhski (pa rallax model). Data taken from 
Ijaroszvnski &: PaczvriskH j2£)02). 



PSPL 




parallax 


2452406.3 ± 0.2 


to 


2452401 


110.8 ± 0.4 




110 


16.751 ± 0.001 


Io 


16.75 


2.370 ± 0.002 


A 


2.249 



good fit to the data using a parallax model, with three more 
fitted parameters, although the result from this more sophis- 
ticated model does not differ greatly from the PSPL model 
presented here (see Table |2J. The anomalous data point is 
around 0.6 mag higher than the fitted PSPL curve, with er- 
rors on this event at ~ 0.01 mag (recorded) and ~ 0.1 mag 
(adjusted). This is a very strong candidate for the detection 
of a planet; the fit suggests a planet with a projected orbital 
radius of around a — I.SRe- A brief survey over q shows 
that this single deviating point is not sufficient to uniquely 
determine a fit - a range of mass ratios between 0.001 and 
0.01 are equally likely. 



3.2.2 2002-BLG-UO 

This event has the largest Ax 2 of all those selected; it is the 
only event to have A\ 2 greater than the A%t = 100 thresh- 
old, roughly corresponding to a 10a deviation. A dark zone 
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Figure 7. Light curve and Ax 2 map for 2002-BLG-055, which is 
the most likely planetary candidate as the strong detection shown 
here is in agreement with independent work. 

on the A\ 2 map (Fig.|SJ, at the image position of the second 
data point after the peak, shows that the q = 1CP 3 binary 
model clearly provides a better fit to the data than a PSPL 
model for this event. A check on the CCD frame responsible 
for the image confirmed that the observation was OK, so the 
point is valid - there are no defects to explain the anomaly 
as non- planetary. However this is a high amplification event 
(Aq = 11.3), and the anomalous point is one of only three 
around the peak. Without a higher sampling rate at this 
critical point, and on the rising part of the light curve, it is 
possible to fit a number of models to the data, so we can 
only regard this event as possible planetary microlensing. A 
preliminary search of a parameter-space involving q, x, and 
y revealed that the highest A\ 2 is found with a mass ratio of 
0.01. The x,y co-ordinates of the 'planet' remained approx- 
imately constant throughout, around the image position of 
the q = 10 -3 anomaly at (1.0,0.58) - corresponding to a 
projected orbital radius of a = 1.16. This is near to Re, as 
expected for a candidate with such a strong influence on the 
event. 



3.3 Lens masses 

Following the method of iTsanras et alJ <|2003T) a lens mass 
can be estimated from the event timescale £e; the best fit 



Figure 8. Light curve and A\ 2 map for the second possible 
planetary candidate, 2002-BLG-140. The data deviates strongly 
from a PSPL model; the peak Ax 2 = 110. This is clearly visible 
as a dark detection zone after the peak. 

mass ratio can give a likely mass for the planetary candi- 
dates. The two events, 2002-BLG-055 and 2002-BLG-140 
have tE = 110.81 and 22.28, corresponding to lens masses 
of 0.53 and 0.24 Mq respectively. The mass ratios then give 
masses of the 'planets' as approximately between 0.6 and 6 
Mj up for 2002-BLG-055, and 2 M Jup for 2002-BLG-140. The 
projected orbital radii, a, for these candidates are 8.6 and 
1.5 AU - if there is a planet around the 2002-BLG-055 lens, 
it is at a distance greater than radial velocity searches can 
currently probe. 



4 THE ABUNDANCE OF COOL PLANETS 

The single good planet-like anomaly in the OGLE-III data 
places a significant constraint on the abundance of plan- 
ets. However, the OGLE-III data are sparsely sampled in 
time, so that many planet-like anomalies can be missed in 
the gaps between data points. We can, however, calculate 
the expected number planet-like anomalies that would be 
detected in this dataset, for various assumptions about the 
planet abundance. 

The total number of detections expected is found by 
summing the individual detection probabilities at each or- 
bital radius - with a in units of Re, so the detection zones do 



Planetary abundance from OGLE 2002 microlensing data. 



7 



Planet Detection for 2002 OGLE III data - 321 events 
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Figure 9. Sum of detection probability, at A\ 2 > 25, 60 and 
100 thresholds, of planets with q = 10 -3 and 10 -4 , as a function 
of orbital radius in Solid lines — q = 10 — 3 , dashed lines - 
g = 10" 4 . 
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Figure 10. Summed probability plot, with a expressed in AU. 
Solid lines — q= 10 -3 , dashed lines - q = 10~ 4 , dot-dashed lines - 
extrapolation to m = Mj up . For each set, the number of expected 
detections is plotted against orbital radius for Ax^ > 25, 60 and 
100. The parabolic shape of the curve is clear; an empirical fit 
fEauation ll7l to it gives an expression for the population of plan- 
ets in terms of a, g, n and A% 2 - dotted line. 



not need to be scaled for lens stars of different masses (Fig. 
|SJ. This was done for mass ratios of 10 -3 and 10 -4 and 
Ax 2 thresholds of Axt = 25, 60, and 100. For isolated data 
points, the detection probability scales with planet mass m, 
so the sum of detection probability for the lower mass ratio 
companions is considerably lower. The peak values of the 
sum of probabilities, and total number of candidates found, 
are displayed in Table [3 

We concentrate primarily on the candidates and detec- 
tion probabilities at Axt = 25, q = 10~ 3 , with a fit includ- 
ing blending. The summed probability peak of ~ 14 (Fig. 

implies that the expected tally of planetary candidates 
should be 14, if every star has a planet with that mass ra- 
tio. The peak probability was found to be at a m Re, as 
expected. 



Table 3. Summary of upper limits on galactic planet population, 
where < n p > is the expected number of planets per star, at a = 
i?E> given n detections where SP were expected, for various values 
of mass ratio q and A\ 2 threshold, Ax|- The figures in brackets 
are the values found if 2002-BLG-140 is taken as a planetary 
candidate. 



<? 


Ax 2 - 


SP 


n 


< ri p > 


10 3 


25 


14.50 


1 (2) 


0.07 (0.14) 


10 3 


60 


10.19 


1 (2) 


0.10 (0.20) 


10 3 


100 


8.15 


(1) 


(0.12) 


10 4 


25 


3.21 


1 (2) 


0.31 (0.62) 


10 4 


60 


2.11 


1 (2) 


0.47 (0.95) 


10 4 


100 


1.63 


(1) 


(0.61) 



< n v > 



The expected number of planets per star is given by: 

(14) 



a 



where n is the number of candidates, and EP represents the 
peak value of the summed probability plots (although this 
result is intuitive, its derivation is non-trivial, and is given 
in Appendix U\l . The q = 10~ 3 analysis gives n = 1, or 
possibly 2 if 2002-BLG-140 is regarded as a candidate, and 
EP = 14.5. This gives < n p >= 0.07, implying that not 
more than 1 in 14 lens stars have planets with q — 10 -3 , 
and a sa IPe ■ 

To look for a higher confidence level of planet detection, 
the analysis was repeated for Axt = 60 and 100. Only 2002- 
BLG-140 has a Ax 2 above the highest level. The predicted 
number of detections is also lower; EP ~ 10 and 8 for Axt = 
60 and 100, respectively. 

Our analysis for q = 10 -4 yielded the same candidates 
with Ax 2 > 25, and a sum of all detection probabilities with 
a peak value of 3.2. This gives an upper limit on the popula- 
tion of less massive planets at 0.31 - implying that smaller 
planets are much more common; 30% of the lens star in this 
sample could have lower mass companions. This confirms 
the expected result from the dN/dm oc m _1 planet mass 
function observed in r esults from radial velocity searches 
jMarcv fc Butlerll2000D . 

Following lTsapr"as et alJ J2003I) . we translate our results 
from dimensionless variables (a/B,E,q) into physical vari- 
ables (a/AU, m/Mjup) using: 



Re = 
and 



1.9AU ( 



V35d/ 



(15) 



(16) 



It should be noted th at Equation ll6l replaces Equation 17 of 
iTsapras et alJ (120031) . which was found to contain an error. 

Fig. 1101 shows that EP(det|a, q, Ax 2 ) is well approxi- 
mated by: 



EP ss 5.5 



io- 



Ax. 

25 



-1/2 



exp ■ 



log (a/4 AU) 



0.4 



(17) 



which captures the scaling of EP with q and (Ax 2 ) 2 , 
and the approximately log-normal distribution of event 
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timescales. Using Equation 1141 this gives an empirical re- 
sult predicting the number of planets per star as a function 
of a, q and A^ 2 , given an observational constraint on n - 
the number of planets detected: 



< n p >^ 3.6x10"°—^ — — exp <^ 3.2 



log! 



I a 



4AU 



.(18) 



Conversion to physical units (Mj up , and AU) broadens 
and lowers the peak in detection capability (Fig. as 
the range in lens star masses gives a range in Re within the 
sampled events. We find a peak in sensitivity at a ~ 4 AU of 
5.5 expected detections; this gives a limit on 'cool Jupiters' 
(q — 10~ 3 ,a ~ 4 AU) of n p ~ n/5.5 = 18%, assuming one 
detection. 



5 COMPARISON WITH PREVIOUS WORK 

The large sample (321 events) studied allows significant up- 
per limits to be placed on the galactic planetary population. 
The result presented here, a limit of 7n percent, is consid- 
erably stronger than t he limits previous yea rs' microlens- 
ing data could provide. iTsaoras et alJ l)2003f) analysed 145 
events, in three years (1998-2000) of OGLE-II data, and 
placed the upper limit at 21n percent, for detection of plan- 
ets with q = 10 -3 at a Axt = 25 level. As the method for 
assessing detection probabilities in the OGLE-II data was 
the same as the one employed by this paper, combination 
of these results was possible; a value of EP « 22 was found 
for all four years of OGLE-II & III (1998-2000, 2002). This 
gives an upper population limit of 5n percent, for the values 
of q and Axt given above. 

Results prese nted by the PLANET collaboration 
jGaudi et al.l2002l) . based on the study of 42 intensely moni- 
tored events, gave a 33 percent limit on Jupiter mass planets 
between 1.5 and 4 AU, at a Axt = 60 level. Looking at our 
result for 'cool Jupiters' at a A\t = 60 level, we find a peak 
at around 4 AU of « 3.5, corresponding to an abundance 
of planets of 28n percent. The results presented are in very 
good agreement with previous limit estimates, improving on 
them due to the larger sample of events studied. These lim- 
its can only be strengthened by inclusion of the more than 
400 events detected in the 2003 OGLE season, which are 
currently under analysis. 



6 SUMMARY 

We analysed the 2002 OGLE-III observations of Galactic 
Bulge microlensing events to assess their capability to detect 
lightcurve anomalies arising from lensing by 'cool' planets. 
From 389 events we select 321 that satisfy criteria for ade- 
quate coverage and compatibility with a point-source point- 
lens light curve model. Our model includes 2 parameters 
that adjust the published error bars to improve their repre- 
sentation of the fit residuals in both faint and bright parts 
of the light curve. We then identify 12 events that appear 
to have significant planet-like anomalies, with a detection 
threshold A\ 2 > 25. Most of these are attributable to iden- 
tifiable causes (e.g. photometry degraded by clouds), leaving 
n = 1 event, 2002-BLG-055 with Ax 2 = 65, which we re- 
gard as a good planetary lensing candidate. 2002-BLG-140, 



with Ax = 110, is also a possible candidate. Our analy- 
sis of detection efficiencies indicates that the 321 selected 
OGLE-III events should reveal 14 planets, with a detection 
criterion Ax 2 > 25, if each of the lens stars has 1 planet with 
a Jupiter- like mass ratio q = 10~ 3 and orbit radius a = Re- 
This implies an abundance of these planets at n v ~ n/14 = 
7n%. Approximate allowance for the range in Re from event 
to event spreads the detection probability over a wider range 
of a, and reduces the peak sensitivity to 5.5 detections for 
planets with a m 4 AU. This gives the abundance of 'cool 
Jupiters' (q = 10~ 3 ,a « 4 AU) per lens star as n p « n/5.5 
= 18n%. The detection probability scales roughly with q 
and (Ax 2 ) _1//2 , and drops off from the peak at a « 4 AU 
like a Gaussian with a dispersion of 0.4 in log a/4 AU (see 
Fie llOl . The peak in sensitivity at a « 4 AU shows that mi- 
crolensing is a complimentary method to other planet search 
techniques, which are more sensitive to close 'hot Jupiters'. 
The anomalies we identify are confined to single data points 
in the OGLE-III data. Higher time-resolution monitoring of 
more events is required to verify the above measurement of 
the abundance of 'cool Jupiters', and to characterise the du- 
ration of the anomalies in order to assess the masses of the 
planets involved. 
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APPENDIX A: PROOF OF UPPER LIMIT 
EQUATION 

The equation for upper limits on planet population 1141 
is found by using Bayes theorem and Poisson statistics 
to describe the probability of detecting planets based on 
a prior assumption about their distribution. The proof is 
given below: 

Definitions: 



Planetary abundance from OGLE 2002 microlensing data. 
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n p = number of planets per star 
E = n p E(P) 

= expected number of detections summed over 
all events 

n = actual number of planets detected over all 
events 



Bayes theorem: 
P(n p \n) = 



P{n\n p )P{n p ) 



(Al) 

(A2) 
(A3) 



Jo°° P (n\n p )P(n p )dn p 

Poisson counting statistics: 

E n e~ E 
P(n\n p ) = — ^ — 

Prior assumption on n p : 
P(n p ) oc n p m , 

m = for prior uniform in n p for n p > 
m = 1 for prior uniform in log (n p ) for n p > 

Expected value of n p given n detections: 

f °° n p P(n p |n)dn p 

< n p \n >= » 

Jo P{n P \n)dn p 

f™n p E n e- E P{n p )dn p 
J °° J B"e- E P(n p )dn p 

Now n p = £/E(P) and dn p = d£/E(P), so 

1 J °° E n - m+1 e- E dE 

< n p \n >= E(P) J~E»-™ e -EdE 

integrating by parts: 

/ x k e~ x dx = -x k e~ x + k x k ~ x e~ x dx 
Jo Jo 

Conveniently, the first term vanishes when evaluated at the 
integration limits x = and x = oo. We therefore have: 



< n p \n > — 



(A4) 
(A5) 

(A6) 
(A7) 



f 

Jo 



E"- m+1 e- E dE = (n - m + 1) / 
/o Jo 

and so: 

(n — m + 1) 



< n„ n > = 



E(P) 



E n-m e -E dE ( Ag ) 



(A9) 



A choice of m = implies an assumed prior on population 
of, for example, 1 planet per star. However this assumption 
gives the same assumed probability of X < n p < X + 1 
for any number of planets per star, X - while lower values 
of X arc intuitively more likely. Taking m = 1 implies an 
assumed probability that is uniform in log(n p ), which has 
the advantage that it automatically gives n p > 0. This gives 
the intuitive result: 

(MO) 



< n p >— 



SP' 



which states that the number of planets per star is given by 
the number found over the expected number. 



APPENDIX B: REJECTED EVENTS 

Events rejected from the analysis as being non-PSPL. The 
rejected events fall into one of four categories, as outlined 



below, depending on which criterion for being an acceptable 
PSPL fit each failed to meet. Remarks on each event from 
the OGLE web page are quoted directly, and marked with 
an asterisk.: 



Rejection categories 



(i) Event listed as a mistake on the OGLE web 
page. 

(ii) Event has over 50 percent error in Ao or £e- 

(iii) Event's oa q and at E are greater than 3 times 
the median. 

(iv) Event has visibly non-PSPL light curve. 



Event 


Category 


Comment 


2002-BLG-002 


(iii) 


Only part of peak observed - parame- 






ters cannot be uniquely determined. 


2002-BLG-003 


(") 


Variable star.* 


2002-BLG-018 


(ii) 


Variable star. Light curve appears 






parabolic in some sections. 


2002-BLG-023 


(ii) 


Variable star?* 


2002-BLG-032 


(iii) 


Only part of peak observed. 


2002-BLG-040 


(iv) 


Variable star?* 


2002-BLG-051 


(ii) 


Double?* 


2002-BLG-057 


(ii) 


Only part of peak observed. 


2002-BLG-058 


(iii) 


Only part of peak observed. 


2002-BLG-059 


(iii) 


Only part of peak observed. 


2002-BLG-068 


(iii) 


Double?* 


2002-BLG-069 


(iv) 


Binary according to PLANET team.* 


2002-BLG-077 


(ii) 


Cataclysmic Variable (I. Bond MOA 






group).* 


2002-BLG-078 


(ii) 


Peak defined by only one data point - 






not enough data to constrain parame- 
ters. 


2002-BLG-080 


(ii) 


Variable star?* 


2002-BLG-081 


(ii) 


Variable star?* 


2002-BLG-089 


(iv) 


Variable star.* 


2002-BLG-090 


(ii) 


Variable star?* 


2002-BLG-099 


(iv) 


Double.* 


2002-BLG-113 


(ii) 


Variable star.* 


2002-BLG-114 


(iv) 


Double?* 


2002-BLG-117 


(ii) 


Second peak well separated from first. 






Wide binary? 


2002-BLG-119 


(iv) 


Double? Variable star?* 


2002-BLG-120 


(ii) 


High Aq event with considerable varia- 






tion in fitted parameters, even on sub- 






sequent runs with same data. 


2002-BLG-127 


(ii) 


Variable star.* 


2002-BLG-128 


(ii) 


Double? Definitely non-PSPL event. 


2002-BLG-129 


(iv) 


Double? Variable star?* 


2002-BLG-135 


(iv) 


Double?* 


2002-BLG-143 


(ii) 


Single? Double?* 


2002-BLG-146 


(ii) 


High Aq event with considerable varia- 






tion in fitted parameters. 


2002-BLG-149 


(iv) 


Severely blended.* Light curve suggests 






lens may be a binary star, as there is a 






visible double peak. 


2002-BLG-151 


(iii) 


Variable star?* 


2002-BLG-152 


(iv) 


Double? Variable star?* 


2002-BLG-159 


(ii) 


Variable star? Peaks throughout data. 


2002-BLG-173 


(ii) 


Very noisy, intrinsically faint, event 






without clear peak. 
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Event 


Category 


Comment 


2002-BLG-175 


(") 


Variable star? Faint, with multiple 






peaks. 


2002-BLG-176 


(iii) 


Variable star? Faint, with multiple 






peaks. 


2002-BLG-194 


(") 


Only part of peak observed. 


2002-BLG-196 


(") 


Parameters vary considerably on dif- 






ferent fits to this small amplification 






event. 


2002-BLG-200 


(") 


Variable star?* 


2002-BLG-202 


(") 


Variable star?* 


2002-BLG-203 


(i) 


Mistake!!!* 


2002-BLG-206 


(ii) 


Very noisy event without clear peak. 


2002-BLG-209 


(") 


Peak defined by only one data point. 


2002-BLG-215 


(") 


Variable star.* 


2002-BLG-228 


(i) 


Artefact from bleeding column* 


2002-BLG-229 


(ii) 


Variable star? No peak visible. 


2002-BLG-255 


(i) 


Mistake!* 


2002-BLG-266 


(ii) 


Any lensing event masked by excessive 






noise. No fit obtainable. 


2002-BLG-272 


(ii) 


Only part of peak observed. 


2002-BLG-273 


(ii) 


Variable star?* 


2002-BLG-278 


(ii) 


No clear peak. 


2002-BLG-284 


(ii) 


High Aq event with considerable varia- 






tion in fitted parameters. 


2002-BLG-286 


(ii) 


High Aq event with considerable varia- 






tion in fitted parameters. 


2002-BLG-307 


(ii) 


Not enough data to define peak. 


2002-BLG-313 


(ii) 


Only part of peak observed. 


2002-BLG-316 


(iv) 


Appears to be either a binary star lens 






or variable star. 


2002-BLG-324 


(ii) 


Only part of peak observed. 


2002-BLG-334 


(iv) 


Strong parallax effect - very long i E , 






event still not complete (2003 season). 


2002-BLG-360 


(iv) 


Only part of peak observed. 


2002-BLG-363 


(ii) 


Only part of peak observed. 


2002-BLG-380 


(ii) 


Only part of peak observed. 


2002-BLG-382 


(iii) 


Only part of peak observed. 


2002-BLG-383 


(iii) 


Only part of peak observed . 


2002-BLG-386 


(ii) 


Only part of peak observed. 


2002-BLG-387 


(iii) 


Only part of peak observed. 


2002-BLG-388 


(ii) 


Only part of peak observed. 


2002-BLG-389 


(ii) 


Only part of peak observed. 



